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Abstract. The peculiar motion of an observer with respect to the CMB rest frame induces 
a deflection in the arrival direction of the observed photons (also known as CMB aberration) 
and a Doppler shift in the measured photon frequencies. As a consequence, aberration 
and Doppler effects induce non trivial correlations between the harmonic coefficients of the 
observed CMB temperature maps. In this paper we investigate whether these correlations 
generate a bias on Non-Gaussianity estimators Jnl- We perform this analysis simulating a 
large number of temperature maps with Planck-like resolution (lmax = 2000) as different 
realizations of the same cosmological fiducial model (WMAP7yr). We then add to these 
maps aberration and Doppler effects employing a modified version of the HEALPix code. 
We finally evaluate a generalization of the Komatsu, Spergel and Wandelt Non-Gaussianity 
estimator for all the simulated maps, both when peculiar velocity effects have been considered 
and when these phenomena have been neglected. Using the value v/c = 1.23 x 10 -3 for our 
peculiar velocity, we found that the aberration/Doppler induced Non-Gaussian signal is at 
most of about half of the cosmic variance a for /jvx both in a full-sky and in a cut-sky 
experimental configuration, for local, equilateral and orthogonal estimators. We conclude 
therefore that when estimating /jvl it is safe to ignore aberration and Doppler effects if the 
primordial map is already Gaussian. More work is necessary however to assess whether a 
map which contains Non-Gaussianity can be significantly distorted by a peculiar velocity. 
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1 Introduction 

The issue of the impact of the peculiar motion of an observer on the Cosmic Microwave 
Background (CMB) is less trivial than it may seem. Its main effect is to induce a large 
dipole, even on a homogeneous map, and actually the value of the dipole is used to measure 
our velocity compared to the CMB rest frame. However the Doppler effect and aberration 
have an effect on all multipoles I when applying a boost to a non-homogeneous map. In 
earlier works [1] and more recently [2-6] this effect has been studied, in order to assess its 
magnitude on two-point correlation functions. In particular in [2-4] it has been pointed out 
that the non-diagonal two-point functions can be used as an alternative method to measure 
our velocity, while in [7] it has been found that the diagonal two-point function (i.e. the 
power spectrum) can be affected leading to a bias that can go up to two standard deviations 
on some cosmological parameters, for a Planck-like resolution experiment. 

It is natural to wonder also if higher-order correlation functions can be significantly 
affected. As stressed already in [1] the aberration and Doppler transformation is a linear 
mixing between multipoles and therefore if the initial spectrum is Gaussian then Gaussianity 
is preserved. However there could still be a sizeable shift on a single random realization 
of a map, especially at large £, since the effect is a mixing of order (5£. More precisely, a 
Non-Gaussianity estimator is a sum over a huge number of combinations of triangles: for 
each triangle there is a large aberration/Doppler induced shift in the measured value of its 
contribution to the total sum, so it is important to check if when summing over many triangles 
the effect becomes negligible for a given experimental sensitivity and resolution. Moreover, 
also the variance of the estimators will be affected by the considered mixing matrix. 

In the present paper we use a direct approach: we simulate 120 maps with Planck- 
like resolution and add to these maps aberration and Doppler effects. We then evaluate 
for each map estimators for the three-point function, both before including aberration and 
Doppler effects and after having accounted for these effects, in order to check how much the 
corresponding Non-Gaussianity parameters can be biased by applying the aberration/Doppler 
transformation on the simulated maps. In particular we focus on the usual /jvx parameters 
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for local, equilateral and orthogonal Non-Gaussianity. We perform our checks both on full-sky 
maps and on cut-sky experimental configurations. 

The paper is organized as follows. In section 2 we shortly review the aberration and 
Doppler effects on a CMB map and on CMB correlation functions. In section 3 we describe 
our procedure to construct simulated maps with and without aberration and Doppler ef- 
fects. In section 4 we define our estimators and show the results. Finally in 5 we draw our 
conclusions. 



2 CMB aberration and Doppler 

The CMB aberration consists in an apparent deflection of the observed CMB photons induced 
by the motion of the observer (e.g. the Planck satellite) relative to the CMB rest frame. The 
difference between the direction n' assigned to the incoming photons by the moving observer 
and the analogous direction n measured in the CMB rest frame is called aberration angle a 
and is give by the relation 

(3sm 2 9 . . 

a ' z =T7^ <"> 

which follows directly from the velocity transformation relating the two frames, assuming 
that the observer moves along the direction z forming an angle 6 with the unit vector n. 
The relative motion between the observer and the CMB also induces a shift in the observed 
frequencies, i.e. the usual Doppler effect. 

Aberration and Doppler induced distortions of a temperature map are quantitatively 
described by the following transformation law 

T'(n') = 7 (1 + /3n • z)T(n) (2.2) 

which relates the temperature T' observed in the moving frame to the temperature T mea- 
sured in the CMB rest frame. In Eq. (2.2) the pre- factor 7(l+/3n-z) accounts for the Doppler 
shift of the photon frequencies and, as a consequence of aberration, the maps T 1 and T are 
evaluated at different directions, namely n' and n which are related by Eq. (2.1). Expanding 
in spherical harmonics both 

T'(n')=^a' lm Y lm (n>) (2.3) 
Im 

and 

r(fi) = £a, ro y, m (A) (2.4) 
Im 

Eq. (2.2) can be rewritten as follows 

j dna M [ 7 (1 + 0h • z)]~ 2 Y$ m ,(n')Y im ,(n) . (2.5) 



/ 



This expression can be then used to compute the desired correlation functions, as we will see 
in the next subsections. 
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2.1 Two-point functions 

Eq. (2.5) can be conveniently rewritten introducing an aberration kernel defined by 

K *t™ = t . n^^ PP^Prf^i:) , (2-6) 

where 



l7 (l-/3x) 1 y ' 1 



and PY l {x) are the associated Legendre polynomials. The harmonic coefficients in Eq. (2.10) 
are hence given by 

^^'<m«l'm. (2.8) 



/ 

l lm 

I 1 



The aberration kernel can be used to write in a compact form all the desired correlation 
functions. Kgi £ m can be estimated: in perturbation theory as explained above, by using 
fitting formulas based on an approximate solution of the integral (2.5) as in Ref. [4], by 
means of the recurrence formula of Ref. [6] or, finally, within the numerical approach used in 
Ref. [7] in the context of cosmological parameter estimation. The aberration kernel allows to 
write the two-point correlation function, namely the covariance matrix, as in the following 
expression 

{ a 'l 1 m a 'l 2 m) = ^ Ki 'i h ™ K \ i 2 m ( a e' 1 m a t' m ) ■ ( 2 - 9 ) 



The integrals in Eq. (2.6) are numerically challenging to compute as a consequence of 
the highly oscillatory character of the integrand which is proportional to the product of two 
spherical harmonics evaluated at different directions. The harmonic coefficients (2.5) can be 
however conveniently estimated by expanding the integrand around f3 equal to zero, being 
our peculiar velocity small in natural units i.e. f3 ~ 1.23 x 1CP 3 . At second order in j3 the 
relation between the primordial ag m in the CMB rest frame and the a' im measured by the 
moving observer takes the following form 

a 'im — (1 + dim)ae m + cj m a^i m + c^ m ae + i rn + dj m an-2m + d^ m a£ + 2 m , (2-10) 
with coefficients cj m , cj m , d £m , dj m , d~ lm given by 

dim = y [(£ + m + 2) G 2 e+hm + l{l - 1) G\ m - l{l + 1) + m 2 - l] 
d L = y K* + 1) Gi+i >m (£ + 2) G i+2 , m ] 

dim = Y \ lG ^ 1 + X ) G ^ra\ (2.H) 

and the function Ge m defined as follows 



£ 2 — m 2 



G ^m = \lj^J- (2-12) 
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We refer the reader to Refs. [1-4] for an explicit derivation of such expressions. From 
Eq. (2.10) one finds for the angular power spectrum measured by the moving observer Ci = 
( a em a tm) following expression 



3(2^ + 1) 



+ s '»-^ 2 imT)}- (2 ' 13) 

Thus, the leading corrections to the primordial angular power spectrum Ci = (ai m a* im ) are 
small, i.e. 0((3 2 l 2 ) (for I <C 1//3), and tend to zero in the limit of flat angular power spectrum, 
when Ci ~ C^+i. This is the reason why in the literature aberration has been neglected in the 
cosmological parameter estimation: it only mildly affects the observed C[ given in Eq. (2.13). 
Eq. (2.13) is however correct only in the full sky approximation, which does not account for 
any foreground contamination of the observed maps. In a more realistic computation, in fact, 
one has to consider that only a fraction of the sky is actually experimentally accessible. This 
can be done by multiplying the original full sky map by an opportune window function to 
mask the background dominated portions of the sky. This procedure leads to a fit of the data 
based on the concept of pseudo-C;, which, as shown in Ref. [7], are significantly more affected 
by CMB aberration and Doppler. Moreover, for t ~ 1//3 the perturbative approach is no 
longer correct. This implies that at large I a more sophisticated description of aberration and 
Doppler effects is required in order to estimate the impact of these phenomena on observable 
quantities, in particular since the Planck satellite very soon will be able to probe the regime 
£ ~ 1//3. In Ref. [7], following a numerical approach similar to the one used in the present 
paper (and explained below), it has been shown that aberration and Doppler can actually 
significantly affect a CMB based cosmological parameter estimation. In that work the impact 
of neglecting Doppler and aberration has been estimated in a Likelihood based analysis of 
temperature maps simulated from a fiducial model including Doppler and aberration. In a 
few simulated maps it was possible to show that neglecting these effects induces a bias on 
the cosmological parameters as large as almost two standard deviations. 

Aberration and Doppler do not only affect the angular power spectrum (see Eq. (2.13)), 
but also the off-diagonal correlation functions, namely the terms (<4 im a^* m )> with £\ ^ 
£2- Notably, the resulting leading corrections are of order f3£ even in a full sky harmonic 
decomposition, contrary to what found instead for the C(_. This property of the off-diagonal 
terms of the covariance matrix has been recently used to forecast the ability of the Planck 
satellite to measure the magnitude and direction of our peculiar motion. In Ref. [3], for 
instance, it has been found perturbatively that a Planck-like experiment will be able to 
measure /3 with an accuracy of about 30%, and the direction of this motion with an error 
of about 20 degree. Ref. [4], instead, using a set of fitting formula for the a^ m , shows that 
such an accuracy would be of about 20% for Planck and even higher for a next generation of 
experiments. 

2.2 Bispectrum 

In Eq. (2.6) we assumed that the harmonic coefficients a' im are measured by an observer 
moving along the direction z coinciding with the z-axis of the harmonic decomposition (2.4). 
Since the direction of our peculiar motion is identified with the direction of the CMB dipole, 
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this choice corresponds to assume that the CMB dipole is aligned along the same z-axis. 
However, when including a mask in the analysis, e.g. to cut the background contaminated 
Galactic plane, it is particularly convenient to work using Galactic coordinates and therefore 
identify the z-axis with the Galactic South Pole-North Pole direction. The latter direction 
does not experimentally coincide with the observed CMB dipole direction which is instead 
characterized by the longitude I ~ 249 and the latitude b ~ 48. This implies that the 
coefficients computed in Eq. (2.5) have to be rotated in order to consistently align the axis 
of the harmonic decomposition with the of Galactic South Pole-North Pole direction. In 
harmonic space this rotation is represented by a Wigner matrix 2^ m /(0, 0, 7), with properly 
chosen Euler angles (ft, 9 and 7. After this Wigner rotation, the harmonic coefficients a» m 
transform as follows 

aim = P w(<M>7)4m' ( 2 ' 14 ) 

-£<m'<£ 

where the coefficients o! im were calculated in Eq. (2.5). The coefficients a£ m are now consis- 
tently expressed in the frame of an observer moving in the direction of the CMB dipole. 

These equations apply to a full sky harmonic decomposition and do not include yet any 
information about the window function W(n) used to account for the partial sky coverage 
of current CMB experiments. If we denote by T(n) a full sky temperature map, its cut sky 
counterpart T(n) takes the following form 

T(n) = W(h)f (n) 

= J2 w em Y em (n)f(h) (2.15) 

I -t<m<l 

where we denoted by W£ m the harmonic coefficients of the window function. Expanding also 
the full sky map 

f(n) = J2^mY lm {n) (2.16) 

lm 

and the cut sky map 

f{n) = Y,ae m Y £m {n) (2.17) 

lm 

one finds 

Q4mi = '/ J '/ J ^mifema %ro 2 (2.18) 
l 2 ~l 2 <m 2 <l2 

where the kernel Ti xmx t iTn% relating the full sky coefficients a^ 2m2 to the cut sky coefficients 
o< imi admits the following representation in terms of the Wigner 3-j symbols 



_ ,„ ( u -m 2 t / (2li + l)(2l 2 + l)(2l 3 + l) 

^imifemj - 2^ 2^ W hm 3 {-L) V 1^ 

tz -h<m s <l 3 



(2.19) 



J \ mi —771-2 W3 

The a,£ m coefficients of Eq. (2.18) were computed in Eq. (2.14). Using now the above equa- 
tions, one can finally derive the three-point function B^^^g^ 3 = (^im%m%m)i namely the 
bispectrum, in the frame of an observer moving along the CMB dipole direction. In terms of 
the kernel J- it reads as follows 

^^Or 13 = X/ XT ^2 ■^imi«m^2m2f'm'^3m3«"tn" (o|m^'m'^"m") • (2.20) 
lm I'm' l"m" 
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3 Simulated temperature maps 



To quantitatively determine the significance of an aberration/Doppler induced Non-Gaussian 
signal, we need first a procedure to effectively simulate a large number of CMB maps observed 
in a reference frame moving along the direction of the CMB dipole. This procedure has been 
developed in Ref. [7] and tested in the context of cosmological parameter estimation. In 
the following we will review this method and its implementation in a modified version of the 
HEALPix 1 code, which we used to calculate the oscillatory angular integrals and the rotations 
in harmonic space introduced in the previous section. These simulated maps will be then 
analyzed in the next section, focusing on possible Non-Gaussian features. As a starting point 

our procedure requires a fiducial cosmological model for which the angular power spectrum 

( f) 

Cy is assumed to be known. The fiducial model assumed in this work is the WMAP 7- years 
best-fit model [10] Given a fiducial model, the map simulation follows five steps relying on 
the following HEALPix routines: 

• synfast: The synfast program generates full sky maps sampling from a fiducial model 

a realization of the ai m coefficients associated with the produced map. These are 

( f) 

drawn from a Gaussian distribution with zero mean and variance equal to the C) . 
We modified this code to allow the user to include the temperature transformation law 
given in Eq. (2.2). We used this modified version of synfast to generate the temperature 
maps studied in the next section. (We restrict ourself to temperature maps, leaving for 
a future work an analysis devoted to aberration effects on CMB polarization maps). 

• anafast: Temperature maps obtained with our modified version of synfast are then an- 
alyzed with the standard HEALPix version of the anafast code. This program expands 
these maps in spherical harmonics calculating the corresponding harmonic coefficients. 
These are the exact (i.e. non perturbative) a' lrn coefficients of Eq. (2.5). 

• alteralm: The a' im coefficients are then rotated in a new reference frame, where the 
peculiar velocity inducing the aberration has the direction of the CMB dipole. To 
perform this rotation we use a modified version of the alteralm code which allows the 
user to choose the Euler angles characterizing the required rotation. The resulting 
coefficients are the ai m given in Eq. (2.14). 

• synfast: A second call to the synfast program generates the temperature map T(n) 
associated with the a^ m coefficients. 

• anafast: Finally, we call again the anafast program to introduce the desired window 
function. 

This procedure provides the CMB maps analyzed in the next section. 
4 Analysis 

An interesting question related to aberration/Doppler effect is whether and how much it 
affects the estimation of primordial non-Gaussianity from inflation. It is well known that 
different inflationary models can produce primordial bispectra with a model dependent shape, 
where the term "shape" is generally adopted to denote the functional dependence of a given 

x http:/ /healpix.jpl. nasa.gov 
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bispectrum on different triangle configurations in Fourier /Harmonic space. So for example 
multi- field models of inflation can produce a "local" shape type of bispectrum, where the 
NG signal in Fourier space is strongly peaked on (ki, &2, £3) triangles with k\ <C &2, ^3 (the 
so called squeezed configurations). Single field models with non-standard kinetic terms in 
the inflaton Lagrangian can instead produce equilateral bispectra, where the signal peaks on 
configurations with k\ ~ k% ~ £3, and so on (see e.g. [16] and references therein for more 
information on bispectra from inflation). For each shape, the strength of the primordial NG 
signal is usually parametrized by a single dimensionless parameter called /nl- A primordial 
bispectrum estimator measures the value of /nl by fitting a given theoretical shape to the 
3-point function extracted from the data. In the idealized case of full sky CMB measurements 
and homogeneous noise, an optimal bispectrum estimator of /nl is derived as (see e.g. [16]) 



^mi,m 2 ,m 3 ^/ NL =l 
h,l2,l:j h,h,ls 



{k,mi} 1 2 3 



&ll ,mi ,rri2 "I3,m3 j 



-.mi,m2,m 3 ^/ NL =l' 

where h lt i 2 j 3 is the theoretical reduced bispectrum for a given inflationary model (obtained 
through a convolution of the primordial shape with CMB radiation transfer functions) and 
^i-ih * s ^ e Gaunt integral. A brute force implementation of formula (4.1) is numerically 
unfeasible for modern satellite CMB experiments (like WMAP or Planck) due to the huge 
number of configurations to evaluate. Komatsu, Spergel and Wandelt showed in [11] that 
when the primordial bispectrum template bi t \ 2 i 3 is separable (i.e. it is a linear combination 
of products of functions of only one of the three I) then equation (4.1) can be conveniently 
factorized and computed in a much smaller number of operations than a brute force evaluatin 
would require. This allow to construct an /nl estimator (from now on KSW) that can be 
applied also to large CMB datasets. Suitable separable approximations have been found for 
the 3 most important shapes in the literature: the local, equilateral and orthogonal shape 
[13]. Following the arguments in e.g. [15] we can then write the KSW estimator for these 
three shapes as: 

/nl P = *Sshape/-^shape j (4-2) 

where -F s hape is the diagonal of the Fisher matrix for the three shapes [11], while S is defined 
as follow: 

Slccal = J r 2 dr J d 2 rv [A(h, r)B 2 (h, r) - 2B(n, r) (A(n,r)B(h,r)) MC 

-A(h,r)(B 2 (h,r)) MC ], (4.3) 

5 equi = -3S local + 6 J r 2 dr J d 2 n{B(n,r)C(n,r)D(ii,r) - B(n,r) (C(n, r)D(n, r)) MC 
-C(n, r) (B(n, r)D(n, r)) MC - D(h, r) (B(n, r)C(n, r)) MC 

-~ [D\n, r) - 3D(n, r) (D 2 (n, r)) MC ] } , (4.4) 
Sortho = 35 cqui - 2 J r 2 dr j d 2 h [D 3 (h,r) - 3D(h,r) (D 2 (n,r)) MC ] , (4.5) 
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and the filtered maps A,B,C,D are denned as: 



A(h,r) = ^fya;(r)(C 1 a)i m Yi m (h) 



(4.6) 



Im 



(4.7) 



Zm 



C(n,r) = ^bmir^C-^mY^n), 



(4.8) 



Zm 



(4.9) 



Zm 



In the formulae above, 6/ is the experimental beam, C is the total power spectrum including 



the CMB signal, Cf MB , and experimental noise, N t , (Q = Cf MB bf + Ni); finally we defined 

a, (3, 7, 5 as: 



bation power spectrum. In the formula 4.3, 4.4 and 4.5, MC denotes Monte-Carlo averages 
over CMB simulations including all experimental features like noise or beam. Those MC 
averages appear in terms that are linear in the a; m 's so that they are zero on average. Such 
terms do not appear in the idealized cubic estimator of formula (4.1). The reason is that 
the formula above defines an optimal /nl estimator only in the case of statistically isotropic 
data. As originally shown in [13], the linear terms have to be included in order to make the 
estimator optimal in the realistic case when statistically anisotropic effects are present in the 
data, in particular sky cut and anisotropic noise. 

In the following we will study the effect of CMB aberration/Doppler on /nl by building 
a KSW estimator and applying it on simulations including aberration/Doppler. The only 
difference between our implementation and the optimal KSW formulae above is that we will 
replace the inverse covariance weighting of the data C~ 1 a with a simple diagonal approxima- 
tion aim/Ci. This will make the estimator only slightly suboptimal when a sky cut is applied, 
while strongly simplifying the numerical implementation. 

4.1 Methodology 

In order to assess the effect of CMB aberration on primordial NG we estimate local, equilat- 
eral and orthogonal /nl from 120 Gaussian maps without aberration/Doppler and compare 
to the same Gaussian maps where the aberration/Doppler is included. This is done in the 
case of full sky and with a mask (KQ75 WMAP7 mask defined in [17]) for an ideal experiment 
(no beam and no noise). 

2 Obtained from a modified version of CAMB code, http://camb.info 




(4.10) 





(4.13) 



gTi is the temperature radiation transfer function 2 and P$ is the primordial curvature pertur- 
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Figure 1. Full sky, local configuration. 



4.2 Full sky results 

Results are shown in figures 1,2 and 3. In the left panel we plot the value of /nl for all the 
shapes and for all the maps. Since the maps are Gaussian the /nl values must be consistent 

with zero within the standard deviations: (7 j*local — 4.2, (7 ^equilateral = 52.3, (7 /.orthogonal = 26.3. 

Jnl /nl Jnl 

That of course holds for both sets of maps (with or without the aberration and Doppler 
effects), as they are both Gaussian. However correlations induced by CMB aberration and 
Doppler can change the estimator variance, and the measured f^L values on a map-by- 
map basis. In blue we plot the /nl values from the Gaussian maps without aberration and 
Doppler effects, and in red the values from the same simulations with aberration and Doppler 
included. The overall effect is small. In order to show the small deviations more clearly we 
plot in the right panels the difference between the /nl obtained from a Gaussian map and 
the /nl from the same map with aberration and Doppler in terms of standard deviations. 
For all the cases the difference is less than half sigma. 

We can conclude that for the case of no mask we didn't find any deviation from the 
case of no aberration and Doppler. 

4.3 Cut sky results 

The same tests performed in section 4.2 are done in this section but with the KQ75 mask 
applied, in order to verify if the presence of a sky cut can enhance the differences. The 
starting Gaussian maps are the same as in the previous section. The presence of a sky cut 

increases the standard deviations to: (7 rlocal — (7 f.cqLii lateral — 67.5, (7 ^.orthogonal = 39.8. 

Jnl J nl Jnl 

Results are shown in figure 4,5 and 6. Even in the cut sky case there is no evidence of a 
significant impact of aberration and Doppler on /nl- 

5 Discussion and conclusions 

In the present paper we have studied the impact of the motion of a peculiar observer, with a 
velocity (3 = v/c = 1.23 x 1CP 3 , on Gaussian maps of the CMB sky with Planck-like resolution 
(lmax = 2000) in order to see if this could significantly bias the measured value of the Non- 
Gaussianity parameter f^L- in fact a peculiar motion induces aberration and Doppler effects 



Figure 2. Full sky, equilateral configuration. 




Map number Map number 



Figure 3. Full sky, orthogonal configuration. 



on the CMB maps, which mixes multipoles with coefficients of order f3£, which are therefore 
0(1) at £ > 800 and so the measured ai m are actually very different from the rest-frame 
ones. Such mixing is still a linear transformation, so on average its effect on 3-point function 
vanishes, but there could be a priori a sizable shift on a single random realization of a map. 
We simulated then 120 maps with Planck- like resolution and we have checked for each map 
how much /jyx can change by applying a boost transformation, using estimators for local, 
equilateral and orthogonal Non-Gaussianity. 

We find that the shift in /jvl is always at most about half standard deviation, both 
for maps with full-sky coverage and with a mask. So, our conclusion is that we can safely 
analyze the map in the boosted frame if our original map is Gaussian, without significantly 
biasing the Jnl measurement. However more work is needed to assess the case in which 
the map is already Non-Gaussian in the CMB rest frame: in such a case there could be a 
considerable shift due to the boost transformation, because the mixing matrix to be applied 
because of aberration and Doppler effects would act on non-zero 3-point functions. 



Figure 4. Cut sky, local configuration. 




Map number Map number 



Figure 5. Cut sky equilateral configuration. 
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